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Abstract 

We study classes of dynamical systems that can be obtained by constructing recursive networks with 
monotone Boolean functions. Stack filters in nonlinear signal processing are special cases of such systems. 
We show an analytical connection between coefficients used to optimize the statistical properties of stack 
filters and their sensitivity, a measure that can be used to characterize the dynamical properties of Boolean 
networks constructed from the corresponding monotone functions. A connection is made between the rank 
selection probabilities (RSPs) and the sensitivity. We also examine the dynamical behavior of monotone 
functions corresponding to filters that are optimal in terms of their noise suppression capability, and find 
that the optimal filters are dynamically chaotic. This contrasts with the optimal information preservation 
properties of critical networks in the case of small perturbations in Boolean networks, and highlights the 
difference between such perturbations and those corresponding to noise. We also consider a generalization 
of Boolean networks that is obtained by utilizing stack filters on continuous-valued states. It can be seen 
that for such networks, the dynamical regime can be changed when binary variables are made continuous. 

1 Introduction 

Biological systems can be viewed as dynamical systems of biomolecular interactions that process information 
from their environment and mount diverse, yet specific responses, striking a balance between stability and 
adaptability. Living systems need to remain stable under variable environmental conditions and yet be able to 
evolve and respond to specific stimuli. This trade-off between adaptability and stability can be captured, in the 
context of dynamical systems analysis, by an order parameter that quantifies the sensitivity of the system to 
perturbations, such as environmental and molecular noise. 

Recent experimental evidence has been mounting in support of a long-standing hypothesis stating that 
living systems operate at the critical regime between ordered and disordered behavior fUHlSlll. Furthermore, 
critical systems have been shown to be optimal in several ways. For example, they are able to store a maximal 
amount of information [5 1 and propagate the information with minimal information loss in a noisy environment 
||6l. Additionally, critical systems exhibit the most complex relationships between their structure and dynamics 
fj]. Thus, optimality of a dynamical system can be considered as a trade-off between noise suppression and 
detail preservation. Such a trade-off is also the central goal of designing digital filters in signal processing. 
Indeed, digital filters can be studied in the context of dynamical systems theory ||8l. 

In this contribution, we consider the well known class of nonlinear digital filters, called stack filters. There 
exists a well established statistical theory for the optimization of stack filters for a given noise distribution 
||9l[T0l. Statistical properties of stack filters can be represented by a set of coefficients and these can be used 
to design an optimal class of stack filters that minimizes some statistical criterion, such as noise variance at the 
output of the filter |[Tn[T2l . These coefficients can also be used to derive statistics that quantify the properties 



of the filters. Since the output of a stack filter is always one of the inputs, rank selection probabilities (RSP) 
quantify robustness and the noise suppression of the filter, specifying the probability that a sample of a particular 
rank in the input window will be used as the output of the filter; while sample selection probabilities (SSP) can 
be used to study the filter's detail preservation, specifying the probability that a particular sample in the input 
window appears at the output |13 |. Thus, these statistics capture the trade-off between robustness and detail 
preservation. For example, the median filter is extremely robust, but preserves details rather poorly, since its 
RSPs are all zero, except the center one, which is equal to one (the one corresponding to the median), while 
the SSPs are all equal (i.e., the uniform distribution). At the other extreme, the identity filter has the highest 
level of detail preservation, but very poor robustness to noise, since the RSPs are all equal, while the SSPs are 
all zero, except for the center one, which is equal to one (the one corresponding to the center position of the 
window). 

We will show an analytical relationship between stack filter coefficients, rank selection properties, and the 
order parameter (essentially, Lyapunov exponent) that characterizes the behavior of the corresponding dynam- 
ical system, thus connecting the statistical optimization of stack filters with dynamical systems theory. 

2 Background 

2.1 Boolean networks and dynamical regimes 

A Boolean network model is a conceptually simple dynamical system model where each node can be in only 
two possible states, on or off. Despite the apparent simplicity, this model class is able to produce highly 
complex behavior, for example, in the form of a phase transition between two dynamical regimes in which 
small perturbations are either attenuated or amplified. A Boolean network model can be defined as follows. 

Let Si{t) € {0, 1}, i = 1, . . . , N, where N is the number of nodes in the network, be the state of ith node 
in the Boolean network at time t. The state of this node at time i + 1 is determined by the states of nodes 
ji,j2,---,jk, at timet as 

S^{t + 1) = fi{Sj,it),Sj,{t),...,Sj,^{t)), (1) 

where /j : {0, 1}'^' — > {0, 1} is a Boolean function of ki variables. A binary vector s (t) = {si (t) , . . . , sn (t)) 
is said to be the state of the network at time t. In the classical model, all nodes are updated synchronously as 
the system transitions from state s (t) to state s {t -\- 1) [Jj. The state transitions are determined by the multi- 
output Boolean function F = {fi, . . ., f^) as s{t + 1) = F{s{t)), where fi is the Boolean function of node i 
with predetermined connections from the nodes ji , j2 j • • • , jki ■ It should be noted that this model can directly be 
generalized to a larger alphabet by defining Si{t) £ {0, . . . , L — 1} and fi : {0, . . . , L — 1}^' — > {0, . . . , L — 1}, 
where L is the size of the alphabet. We restrict our attention to the case L = 2. 

To construct a Boolean network, the inputs ji , , • • • > jki for each node i needs to be determined. This can 
be done by selecting the inputs randomly among all N nodes or by selecting the inputs using some systematic 
pattern. The number of inputs ki can be endowed with a probability distribution, such as the power-law |[T4llT5l 
or Poisson distribution, with a mean K = E [ki] , known as the average connectivity of the network. 

Once the connections have been set, we can choose a Boolean function fi for each node. Functions can be 
parameterized by the bias h = E [fi], the probability that the function outputs one on an arbitrary input vector. 
If 6 = 0.5, then the function is said to be unbiased. The functions can be selected randomly among all 2 ' 
Boolean functions or they can be selected from some class of functions |[T6l[T7l[T8l[T9l . If both the functions 
and connections are selected randomly, then the obtained network is called a random Boolean network (RBN) 

m. 

Since a Boolean network is a discrete system, it has a finite state space. Thus, every state trajectory, that 
is, a path through the state space from any initial state, will eventually return to one of the previously visited 
states. This kind of a state cycle where the same states are repeated infinitely is known as an attractor cycle and 



the states within the cycle are called attractor states. A set of states that leads to the same attractor is called the 
basin of attraction ll20l . 

Boolean networks, as models of dynamical systems, can operate in the ordered or chaotic regimes, or at 
the phase ti^ansition boundary between these two regimes [21 ]. This phase transition regime has been referred 
to as the edge of chaos [1]. When a network is operating in the ordered regime, it is intrinsically robust while 
its dynamical behavior is simple. The robustness can be observed through both the structural and transient 
perturbations. Perturbations have a small effect on the behavior of the network. Networks in the chaotic 
regime, on the other hand, are extremely sensitive to perturbations. Even a small perturbation will quickly 
propagate through the entire network. Thus, networks in the chaotic regime are not robust in that they are 
not able to coordinate macroscopic behavior under perturbations. A phase transition between the ordered and 
chaotic regimes represents a trade-off between the need for stability and the need to have a wide range of 
dynamical behavior to respond to variable perturbations [l]. 

By varying the parameters K and b in the random Boolean network model, a dynamical phase transition 
can take place. The average sensitivity of the Boolean network 

a = 26(1 - b)K (2) 

can be used to determine the dynamical regime. If it > 1 then the system is chaotic and for u < 1 the system is 
ordered Il22ll23ll24ll25i . It is easy to see that for unbiased random Boolean networks, the critical connectivity 
is Kc = 2. 

In a sense, critical systems are maximally responsive to the useful information in their environment while 
being able to reliably execute their behaviors in the presence of uninformative variation in this environment. A 
hallmark of critical behavior is the spontaneous emergence of complex and coordinated macroscopic behavior 
in the form of long-range spatial or temporal correlations. Such coordination across many scales enables 
information to propagate over time from one part of the system to another with a high degree of specificity 
and sensitivity. These aspects of criticality support the idea that living cells, as complex dynamical systems of 
interacting biomolecules, are dynamically critical systems. 

2.2 Stack Filters 

Stack filters constitute an important class of nonlinear filters based on monotone Boolean functions |[26l . Statis- 
tical properties of stack filters have been studied in terms of output distributions and moments for independent 
and identically distributed input signals lITlfTTI . Consequently, it becomes possible to optimize stack filters in 
the mean square sense. In other words, the knowledge of the input distribution allows one to find a stack filter 
or a set of stack filters that minimize the output variance. 

Rank selection probabilities (RSP) and sample selection probabilities (SSP) are probabilities that the output 
equals a sample with a certain rank and certain time-index in the filter window, respectively. The output 
distribution of a stack filter can be expressed in terms of its RSP's. On the other hand, SSP's give us information 
about the temporal behavior of stack filters. This information is important for examining the detail preservation 
properties of stack filters. 

Let a = (ai, • • • , a„) and /? = ■ ■ ■ , (3n) be two different n-element binary vectors. We say that a 
precedes /?, denoted as a ^ /?, if Oj < /3j for every i, I < i < n. If a -/^ (3 and j3 ^ a, then a and j3 are 
said to be incomparable. Relative to the predicate -<, the set of all binary vectors of a given length is a partially 
ordered set. A Boolean function / {xi, ■ ■ ■ ,Xn) called monotone if for any two vectors a and /3 such that 
a ^ /3, we have / (a) < /(/?). The class of monotone Boolean functions is one of the most widely used and 
studied classes of Boolean functions. 

Let i?" denote the Boolean n-cube, that is, a graph with 2" vertices each of which is labeled by an n- 
element binary vector. Two vertices a = (ai, • • • , a„) and /3 = {Pi, - ■ ■ , are connected by an edge if and 
only if the Hamming distance p{a, (3) = Y^^=i {'^i ® A) = 1) where © is addition modulo 2 (exclusive OR). 



The set of those vectors from E"' in which there are exactly k units, < A; < n, is called the kth layer of i?" 
and is denoted by E"-'^. The Hamming weight of the vector x G {0, 1}" is w{x) = p{x, 0), where is the 
all- zero vector. 

A stack filter is defined by a monotone Boolean function. A continuous stack filter Sf{-), based on function 
/, is obtained by replacing conjunction (•) and disjunction (+) operations by min(-) and max(-) operations, 
respectively. For example, the monotone Boolean function f{xi,X2,x^) = X2 + xix^ corresponds to the 
continuous stack filter Sf{Xi, X2, X^) = max(X2, min(Xi, X3)). 

Suppose that the input variables of some stack filter Sf{-) are i.i.d. random variables with distribution 
F{t) = Pr{Xi < t}. Then, it is well known |27| that the output Y = Sf{Xi, . . ., X„) of the stack filter has 
output distribution 

n-l 

a>(t) = ^A,(i-F(t)rF(tr- (3) 

i=0 

where 

A = \{xeE^'':f{x) = 0}\. (4) 

It is known 1^281 |29l that rank selection probabilities = Pt{Y = where is the ith order 

statistic, are obtained as. 

Pi ~ ( n \ ~ 1 n Y' ^ 

\n—i) \ri— 1+1/ 

It can also be shown that the output distribution function of the filter can be expressed in terms of the RSPs as 

n 

$(t) = ^p,F(,)(t), (6) 

i=l 

where is the cumulative distribution function of the ith-order statistic for i.i.d. inputs 111. 

Let Sj = Pr{y = Xj} be the jth SSP and s = (si, . . ., s„) be the vector of SSPs. In |l30]|, it was shown 

that Sj = dj{l) — dj{0), where 



(7) 



and f-\k) = {xeE^ : f{x) = k}. 

For a given input distribution, the mean square optimization of stack filters can be performed as follows. 
The variance fj,2 = E{(Y — E{Y})'^} of the output Y of the stack filter can be written as 

n-l /n-l \ 2 

^i2 = Yl ^^^(^' 2, n, i) - AiM{F, 1, n,i)] , (8) 

i=0 \i=0 / 

where 

r°° . d 

(9) 



/oo J 
x^^{{l- F{t)yF{t)''-')dx. 



Thus, for a given noise distribution, the goal of optimization is to find parameters Ai (equivalently, the rank 
selection probabilities) such that the objective function ([8]) is minimized. An optimal stack filter, specified by 
monotone Boolean function /, is then constructed using the obtained parameters A^. As a result of optimization, 
we obtain a set of parameters Ai that define a class of stack filters. In a given class, all stack filters are 
statistically equivalent since they all possess the same parameters A^. Thus, they will also have the same rank 
selection probabilities p,. An approach to select the best filter, in terms of its ability to preserve details, from 
the class of such statistically equivalent optimal stack filters has been proposed in ||T3t . 



2.3 Average sensitivity of Boolean functions 

Let / : {0, 1}" {0, 1} be a Boolean function of n variables xi, . . . , Xn- Let 

df (x) /dxj = f (x^^^^A © / (x^^^^A (10) 



be the partial derivative of / with respect to xj, where x^^'''^ = (xi, . . . , xj^i, k, Xj+i, . . . A; = 0, 1. 
Clearly, the partial derivative is a Boolean function itself that specifies whether a change in the jth input causes 
a change in the original function /. Now, the activity of variable Xj in function / can be defined as 



«^=4 E 9fix)/dx,. (11) 



Note that although the vector x consists of n components (variables), the jth variable is fictitious in df (x) / dxj . 
A variable Xj is fictitious in / if / (^x^^'^^) = f [x^^'^^) for all x^^''^^ and x^^'^\ For an n-variable Boolean func- 
tion /, we can form its activity vector = a(, . . . ,an . It is easy to see that < < 1, for any 

j = 1, . . . ,n. In fact, we can consider aj to be a probability that toggling the jth input bit changes the 
function value, when the input vectors x are distributed uniformly over {0, 1}". Since we're in the binary 
setting, the activity is also the expectation of the partial derivative with respect to the uniform distribution: 

aj = E[df {x)/dxj\. 

Another important quantity is the sensitivity of a Boolean function /, which measures how sensitive the 
output of the function is to changes in the inputs. The sensitivity {x) of / on vector x is defined as the number 
of Hamming neighbors of x on which the function value is different than on x (two vectors are Hamming 
neighbors if they differ in only one component). That is, 

{x) = |{i e {1, . . . , n} : / (x © ei) / / (x)}| (12) 

n 

= ^x[f{x(Be,)^f{x)], 
1=1 

where is the unit vector with 1 in the ith position and Os everywhere else and x [^] is an indicator function 
that is equal to 1 if and only if A is true. The average sensitivity is defined by taking the expectation of 
S'^ (x) with respect to the distribution of x. It is easy to see that under the uniform distribution, the average 
sensitivity is equal to the sum of the activities: 



It 

E [sf (x)] = Y^E[x[f{x®e;)^f (x)]] (13) 

n 



1=1 

i=l 



Therefore, s-^ is a number between and n. For a random Boolean network, the mean of the average sensitivities 
S'^* of the Boolean functions fi, i = 1, . . . ,N, corresponding to each of the iV nodes is the network average 
sensitivity a given in (|2]| |[25l . 



3 Average sensitivity and rank selection probabilities 

As defined in (fT3] ). to compute the sensitivity of a monotone function / with n variables, we need to go through 
all the one-bit changes (from each input vector x G -E") and see how many times, doing this, the output 



changes. We can divide the task into n sub-tasks, in each of which we only consider bit changes (perturbations) 
from a vector x with w{x) = k into a vector y with w{y) = k + 1 (k = 0,1, . . .,n — 1). We can denote this set 
of perturbations {x, y) by Fk for simplicity. 

For each of the (^) vectors x G E"-'^, it holds that 

\z : {x, z) G Fk\ = n — k, 

and conversely, for each vector y with w{y) = k + lit holds that 

\z : iz,y) £Fk\ = k + l. 

In total. 



IF. 



(n — k) 



n 
k + 1 



(k + l), 



(14) 



so that there are this many perturbations to check in this sub-task. Since we know the total number F^ = 
n X 2"~^, we can just as well look for bit changes in each F^ that do not change the output of /. 

Based on the monotonicity of /, we know that for each x G E'"-'^ and /(x) = 1, all perturbations {x, y) in 
Fk will not alter the output of the function (since f{y) = l as well). There are — such states and hence 



A, 



(n — k) 



(15) 



such non-output-altering perturbations in F^. 

On the other hand, we know that for each y G ]^n,k+i j^^-j _ perturbations (x, y) will also be 
non-output-altering (since f{x) = due to monotonicity). There are A^j^i such states and hence 



Ak+i{k + l) 



(16) 



such perturbations in F^. 

The number of output-altering perturbations in F^ can now be computed by taking out the contributions of 
([T5] ) and ([T6l ) from the total number in (fT4l) . i.e. 



n 
k + l 



{k + l) 



A, 



{n-k)- Ak+i{k + l). 



Summing together the contributions for all k and normalizing with 2"~^ (the total number of vectors from 
which a perturbation may be made, divided by 2 to take out symmetric perturbations that have not been counted 
separately) gives the sensitivity of / as 



n-l 



n 



—y 

n-l 

k=0 
^71—1 /n 



{k + l)Ak+i + 



Ak]{n- k) 



(17) 



Taking into account the fact that X]fc=o (fc) {n — k) = n x2'^ ^ , this can be simplified to obtain 



n-l 



^ = ^ E [^-^(^ -k)-{k + l)Au+i] , 



k=0 



which shows the relationship between the coefficients A^ and sensitivity S'^. 
From we can solve 



[Akin -k)- Ak+i{k + 1)] = pn-k [^]{n-k). 



n 



Inserting this into dTT] ) results in 



n 



n-l 

tE 

fc=0 



n 



(n - k) -Pn-k{ , ]{n-k) 



^ n— 1 
n-l ^ 



fc=0 



1 " 

n-l 



fc=l 



n 



n 



A: 



giving the connection between the rank selection probabilities pi and sensitivity s-^ . 

Interestingly, it has been shown that almost all stack filters are robust in the sense that all but the cen- 
tral four rank selection probabilities are nonzero |[29l . Moreover, the central two rank selection probabilities 
asymptotically dominate over the outer two rank selection probabilities. Using similar derivations, it is possi- 
ble to compute the average sensitivity of a so-called typical monotone Boolean function. The size of the set of 
typical monotone Boolean functions asymptotically approaches the total number of monotone Boolean func- 
tions. Equivalently, the probability that a randomly picked monotone Boolean function is a typical function 
approaches unity. Because of the super-exponential growth of the number of monotone Boolean functions, the 
asymptotic convergence occurs very rapidly with excellent accuracy even for as few as seven input variables. It 
can be shown that the expected average sensitivity of a typical monotone Boolean function is given by 



n2~ 



n 



n/2 



^ (2-"/2-i + 1 



(18) 



for even n.The case of odd n is similar to derive, but its expression is much more cumbersome (see f3Ti)- The 
expected average sensitivity of a typical monotone Boolean function is plotted in Figure [T] It can be seen that 
stack filters with as few as five inputs are already chaotic, in terms of the average sensitivity being interpreted 
as a dynamical order parameter 



4 Stack filter performance in terms of dynamical behavior 

To relate the noise suppression properties of a stack filter to the dynamical behavior, captured by the average 
sensitivity, we filtered various test signals and noise distributions with all five-input stack filters. Figures |2] 
and [3] show the mean square error (MSE) as a function of filter sensitivity for blocks and Heaviside signals, 
respectively. In addition. Figured shows the MSE performance for different degrees of salt and pepper noise. 

Based on these results, it is clear that filters that are chaotic, in terms of their sensitivity, perform best 
under various noise conditions. Simulations with salt and pepper noise show that as more noise is added, the 
performance of ordered and critical filters deteriorate more that the chaotic ones. 

In light of what we known about the dynamics of random Boolean networks, this observation is somewhat 
surprising. A large body of work with Boolean networks has established that dynamically critical systems are 
the ones that strike the optimal compromise between robustness to noise and reliable propagation of informa- 
tion. Thus, it would be expected that dynamically critical stack filters would have performed optimally under 
variable noise conditions in terms of noise suppression and detail preservation. 

There are several aspects of signal denoising that are essentially distinct from the analysis of dynamical 
systems. The theory, which relates the average sensitivity of a Boolean network to the dynamical regimes, 
is based on the propagation of small perturbations. That is, to quantify the dynamical behavior, we measure 
whether small perturbations that are introduced into the state of the system are amplified or attenuated, leading 
to chaotic or ordered behavior, respectively. 

Noise filtering is fundamentally different in terms of what we know about propagation of small perturba- 
tions. A noisy signal is frequently the result of a large perturbation that has affected the value of majority of 




Figure 1: The expected average sensitivity of a typical monotone Boolean function versus the number of 
variables n. The cases of even n (solid line) and odd n (line with circles) are shown separately. 




Figure 2: Mean square error versus sensitivity of stack filter. Data is shown for blocks signal under (a) Gaussian 
and (b) bimodal noise (see fTSl for details on noise distributions). 




Figure 3: Mean square error versus sensitivity of stack filter. Data is shown for Heaviside signal under (a) 
Gaussian and (b) bimodal noise. 
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Figure 4: Mean square error versus sensitivity of stack filter. Data is shown for Heaviside signal under salt and 
pepper noise with probability (a) 0.001, (b) 0.01 and (c) 0.1. 
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Figure 5: Information theoretical order parameter (Derrida curve) computed for continuous signals, (a), (b) and 
(c) show filter classes that are ordered, critical and chaotic in terms of the average sensitivity of the Boolean 
function implementing the stack filter, respectively. All K = 5 stack filters are shown. D{t) is the information 
distance between two input signals and D{t + I) is the distance between filtered signals. Each line is the 
averaged behavior of a single stack filter. See ||4J for the details on the construction of the Derrida curve. 

the signal values (consider additive or multiplicative noise that affects every signal value or pixel). Thus, the 
insights derived from the small perturbation analysis of dynamical systems do not necessary hold. The mean 
squared error of the filter does not necessary capture the trade-off between noise suppression and detail preser- 
vation that is observed when small perturbations are studied. Very little is known about the dynamical systems 
behavior when only the observations about the response to large perturbations are available. 

4.1 Stack filter as a generalization of a monotone Boolean network 

An important aspect of stack filters is that they can be viewed as a generalization of the monotone Boolean 
network model to a continuos model. As a stack filter corresponds to a monotone Boolean function and filtering 
can be implemented in the form of a Boolean network, the stack filter forms a network that can be run for 
continuos signals. This insight gives an interesting possibility to analyze the relationships between Boolean 
and continuos models. 

Recently, there have been several attempts to link the properties of Boolean networks to the behavior of 
classes of continuous or more detailed discrete models. For example, the existence of dynamical regimes has 
been demonstrated in models other than Boolean networks. The interpretation of stack filters as generalizations 
of Boolean networks allows us to study the connections between Boolean and continuos models directly. 

In recent work (4], we introduced an information theory based order parameter that can, in principle, be used 
to quantify the dynamical behavior of any model class, discrete or continuous. This parameter measures the 
average tendency of the dynamical system to attenuate or amplify informational distances (computed using real 
world compressors) between states of the system. We utilize this order parameter to quantify the dynamical 
behavior of stack filters in both the binary and continuos forms. In the binary case, we introduce single bit 
perturbations to states of the system with varying probability. In this standard analysis, the information theory- 
based order parameter has been shown to accurately capture the dynamical behavior, characterized by the 
average sensitivity of the Boolean function 0. For the continuous case, we can introduce noise in the form of 
salt and pepper noise with varying probability. Results for the continuous case are shown in Figure [5] 

Remarkably, this analysis suggests that the dynamical regimes of the networks are different when the input 
data are changed from binary to continuous. For example, a Boolean network that is chaotic in terms of 
average sensitivity (and confirmed with the information distance -based order parameter) may not be chaotic 
in the continuous domain in terms of the same information distance -based order parameter. This suggests that 
the properties that are derived for the monotone Boolean functions do not necessary hold when the function is 



generalized to the continuous case. 



5 Conclusions 

We have shown an analytical relationship between stack filter design coefficients Ai and the average sensitivity 
of the monotone Boolean function. In addition, we presented a formula for the average sensitivity of typical 
monotone Boolean functions. This formula implies that typical stack filters are dynamically chaotic. 

These results suggest a novel perspective that can potentially be utilized in the design of stack filters for 
specific filtering tasks. It ties the design of stack filters to a more general dynamical systems framework. Thus, 
we can also utilize stack filter optimization algorithms to design Boolean networks that are statistically optimal 
under given noise distributions. It will be of interest to compare the properties of such statistically optimal 
ensembles of networks to the properties of other ensembles under different noise conditions. By estimating the 
noise distributions from biological data and designing optimal Boolean networks may also help us gain insight 
into the dynamical behavior of biological systems. 

We also studied the behavior of stack filters under various noise distributions by relating the MSE to the 
average sensitivity of the filter. This analysis suggests that chaotic filters perform best under various noise 
conditions. This is contradictory to what was expected based on our knowledge about Boolean networks. 
However, subsequent analysis shed Ught on this observation. Dynamics of Boolean networks are defined under 
the assumption of small perturbations while the filtering of noisy signals is a fundamentally different problem. 
In addition, our parallel analysis of stack filter dynamics in Boolean and continuous cases suggests that the 
dynamical behavior of monotone Boolean functions does not directly generahze to continuous systems. 

Future work should focus on understanding these aspects in more detail. More work needs to be done 
to be able to quantify the dynamics of a system from large perturbations. It will also be of interest to study 
the connections between dynamical behavior of Boolean and continuous models. Stack filters will be a useful 
model class for these studies due to the availabihty of direct generalizations from the Boolean to the continuous 
domain. 
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